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Abstract 

It is shown that generalized CDT, the two-dimensional theory of quantum gravity, 
constructed as a scaling limit from so-called causal dynamical triangulations, can be 
obtained from a cubic matrix model. It involves taking a new scaling limit of matrix 
models, which is more natural from a classical point of view. 

Introduction 

The great versatility of matrix models or matrix integrals in theoretical physics is 
well illustrated by their particularly beautiful application in two-dimensional Euclidean 
quantum gravity (see [1, 2, 3, 4] for reviews). This theory can be defined as a suitable 
sum over triangulations, so-called "dynamical triangulations" (DT), whose continuum 
limit is obtained by taking the side lengths a of the triangles to zero. The method 
of DT was originally introduced as a nonperturbative worldsheet regularization of the 
Polyakov bosonic string [5, 6, 7]. There it was used with success (or to disappointment, 
depending on ones taste) to show rigorously that a tachyon-free version of Polyakov's 
bosonic string theory does not exist in target space dimensions d > 1 [8]. However, 
when viewed as a theory of 2d quantum gravity coupled to matter with central charge 
c < 1, the theory - noncritical string theory - is perfectly consistent, and matrix mod- 
els have been used to solve in an elegant way the combinatorial aspects of the DT 
construction, where one sums over random surfaces glued together from equilateral 
triangles. 

The DT approach possesses a well-defined cut-off, the length a of the lattice links. 
As has been discussed in many reviews (for instance the ones mentioned above), a con- 
tinuum limit can be defined when the lattice spacing is taken to zero while simultane- 
ously renormalizing the bare cosmological constant and possibly other matter coupling 
constants. However, the continuum limit in question has some unconventional proper- 
ties. We will show that there is another way of taking the scaling limit of the matrix 
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Pi 

Figure 1 : The gluing of triangles via Gaussian integration. 

models, which still relates them to a summation over triangulated random surfaces, 
the so-called causal dynamical triangulations (CDT) [9]. We will show that this limit 
is in a way more natural and it corresponds to starting out from a "classical" matrix- 
model theory. In accordance with this it does not lead to the somewhat unconventional 
renormalization encountered in the standard DT approach. 



Hermitian Matrix Models 

We can define the Hermitian matrix model for N xN matrices as a formal power series 
in 9 



Z(g) = /d^e-^^M* 3 ) (1) 
dip = TT d~Re(f> a p Y[ dIm0 Q( g. (2) 



a</3 a</3 



The integral can be evaluated in the standard way by performing all possible Wick 
contractions of (Tr 4> 3 ) k and using 

(<t>ap<t>a>p>) =C / # e - * lMa <t>c0<l> a 'P' = <V<V, (3) 



The geometric interpretation in the context of DT is illustrated in Fig. 1 . An index 
is assigned with each vertex in a triangle and a matrix <f> a p is assigned to the link (the 
side in the triangle) which contains the vertices labeled a and j3. In that way we can 
associate tr 3 with a triangle and the term (</> 3 ) fc in (1) with k triangles. Performing 
the Gaussian integrations indicated in eq. (3) we are gluing together all triangles in all 
possible ways by identifying links as illustrated in the figure. This way of performing 
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the integral will result in an asymptotic power series in g which can be Borel summed 
for g < 0. 

If we want to perform the integral without power expanding the c/> 3 part of the 
potential one takes advantage of the invariance of the action under <j> — > U 4>W , where 
U € U (N), Thus the action depends only on the eigenvalues li of cf> and we can make 
the following decomposition of the measure of integration: 

d<j> e-i^W oc dU(N) f[ d£i (T^U) [] % - (4) 

i=l i<j 

where Yi%<j ~ ^j\ 2 * s tne Jacobian, changing from <fi to its eigenvalues and the 
unitary matrix U. The integral over the U matrices is now trivial since the action is 
independent of U, 

The "classical" limit is obtained for g s -4 0, where all eigenvalues are lumped 
together at £o, where 

V'(£ ) = 0. (5) 

However, for g s > the integration over the non-diagonal matrix elements pro- 
duces the Vandermonde determinant Yii<j ~ ^il 2 ' which acts as a "quantum" cor- 
rection, a repulsion between different eigenvalues. The result is that eigenvalues are no 
longer concentrated at £q. In the large N limit they occupy an interval around 1$. 

As an example we let us consider the matrix potential 

1 t r f js\ 1 ( „x , 1 J.2 9 ,3 



g s g s \ 2 3 / 

shown in Fig. 2. It is clear that this specific matrix integral only exists as the formal 
power series for finite N. For infinite N the eigenvalues will condense in an interval 
around £q. This interval is determined by a large N saddelpoint equation. We will not 
here discuss the solution to that equation, but only mention that the so-called resolvent, 
which determines the distribution of eigenvalues is given by: 

/I. 1 \ !/",,! 1 -JV trV(4>) 



For V(4>) = —g<fi + \§ — |0 3 one has from the saddelpoint equation or by other 
methods, like the so-called loop equations, to leading order in N: 

w(z) = i (V'(z) + g(z - bW(z - c)(z - dfj , (8) 

where the constants b, c and d are determined by the requirement that w(z) — > 1/z 
for z — > oo. It should further be noticed that in the large N expansion each term has 
analyticity structure like w(z) in the complex z-plane, i.e. a branch cut between c and 
d. 
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Figure 2: The graph — g<j> + \<\> 2 — |</> 3 



The conventional scaling limit 

The usual scaling limit of the matrix model, relevant for non-critical strings and 2d 
Euclidean quantum gravity coupled to matter, is obtained (for fixed g s ) by fine-tuning 
g such that b(g) = c(g). At this point the analytic structure of w(z) changes from 
(z — c(g)) 1 / 2 — > (z — c(g)) 3 / 2 , and this change can only be accommodated by invoking 
arbitrary high k in the sum 

This is why one geometrically can imagine a "continuum" limit where the size of each 
triangle shrinks to zero while the continuum size of the surface stays constant. More 
precisely, for our specific model one has: 

g = g c (l - Aa 2 ), z = c(g c ) + aZ, a -> 0, (10) 

w(z) = I (y'(z) + gVc(g c ) - d(g c ) a 3 / 2 W E (Z, A) + 0(a 5 / 2 )) , (11) 
where the "continuum disk amplitude is 

We(Z, A) = (Z - v/2A/3) yf Z + 2^2k/Z. (12) 

Here a has the interpretation as the length of the side of the triangles (polygons) which 
appear in V(^). 

Notice that actually the non-scaling part V'(z)/2 dominates when a — > and ren- 
ders the average number of polygons present in the ensemble with partition function 
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w(z) finite, even at the critical point. This somewhat embarrassing fact can be circum- 
vented by differentiating w(z) a sufficient number of times with respect to g and z, 
after which these "non-universal" contributions vanish since they are polynomials in g 
and z, but for the disk amplitude itself there is no such escape. 



The new scaling limit 

Until now we have considered two limits: the "classical limit": g s = and "conven- 
tional scaling limit" of non-critical string theory: g s > and g — > g c (gs)- Is it possible 
to find a new, non-trivial scaling limit, closer to the classical limit when g s — > 0. The 
answer is yes [10]. 

If one works out the details close to the conventional critical point b(g c ) = c(g c ) 
one has 

3c(.9 s ) = i(l-|9 s 2/3 + 0(.g^)), (13) 

zM = c(g c , g a ) = l + gl /z + 0( 5s 2/3 ), (14) 
c(g c )-d(g c )=Ag]l 3 + 0{g 2 J 3 ) (15) 
A non-trivial scaling can now be obtained if we fine-tune g s — > as 

g s =G s a 3 . (16) 

Again the scaling parameter a can be given the geometric interpretation as the link 
lengths of the polygons in V(4>). Note that the length of the cut goes to zero as a — >• 0, 
thus we are closer to the "classical" limit. However, it will survive in the continuum 
limit: 

g = g c {g s )(l - a 2 A) = g(l - a 2 A cdt + 0(a 4 )) (17) 
z = z c + aZ = z + aZ c( it + 0(a 2 ) (18) 

A cdt = A+^G 2 J 3 , ff=i Z cdt = Z + Gl/ 3 , z = l. (19) 
Using these definitions one computes in the limit a — > that 



I A cdt - \Z 2 cdt + \{Z cAt - H)^J{Z cdt + Hf 



4G S 
H 



w(z) = * . (20) 



(2A cdt )3/2 



0, h = H/^2A cdt (21) 



w (z) = - W cdt (Z cdt , Adt, G„) (22) 
a 

Thus we have a situation where, contrary the situation we encountered taking the con- 
ventional scaling limit discussed above, the disk amplitude w(z) really scales as one 
would expect from an ordinary correlation function. 
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We can now take the limit G s — > and we obtain 



^ cdt (Z cdt ,A cdt ,G s ) ->■ (23) 

^cdt + V ^ A cdt 

Thus we see that the cut where the eigenvalues are located shrinks to a point, indicating 
we have precisely the classical situation discussed above. This is indeed the case as we 
will discuss further shortly. 

If we on the other hand take the limit G s — > oo we obtain for the square root part 

of W c dt(^cdt, A cdt ,G s ) 



(Z cdt -H)J(Z cdt +H) 2 -^ 



2G S 



G s ~ 5/6 (Z-V2A/3) ^Z+2y/2Aji 



Thus we recover the standard continuum disk function We(Z,A) from (12) times a 
factor GJ 5 ^ 6 . If we write g a = a 3 G s (in accordance with (16)) and keep g s constant, 
while taking a — > 0, it means that G s — > 00 as a~ 3 . Thus G s ~ a 5 / 2 and we 
precisely recover the square root term in (11) if we remember that w(z) and Wcdt 
according to (22) differ by a factor a. However, the part not related to the square root 
in W c dt will not scale in the limit G s — s- 00, in accordance with the previous discussion 
of standard scaling related to (1 1). 



The matrix model 

The new scaling can be obtained by a simple change of variables: 

<f)^ z I + a<5> + 0(a 2 ) (24) 
Up to a independent term we then have 

A n$> — i$ 3 

vw = y(*), n*)= 2Gs 6 , (25) 

where V(4>) is the potential given by (6). Thus we can write 

Z(g,g s ) = a N2 Z(A cdt ,G s ), Z(A cdt ,G s ):= / d<t> e~ N ^(*) 



(26) 



The change of variable (24) explains immediately and in a trivial way the scaling (22) 
of w(z): 

— = - „ 1 . w(z) = -W cdt (Z cdt ,A cdt ,G s ) (27) 

z- <p a Z cdt — <P a 

This relation is obviously correct to all orders in N What is truly surprising is that the 
new scaling limit is itself a matrix model defined by V($>) [11]. The continuum limit 
a — > is thus described by a matrix model. This bear some resemblance with the 
Kontsevich matrix model, even the matrix potential is somewhat similar. But contrary 
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Figure 3: The potential V(£), (25), with the local minimum £q — 



— \/2A c dt- 



to that model, where the continuum objects are the modular spaces of Riemann surfaces 
it turns out that the present matrix model actually describes a set of "real" continuum 
surfaces as we will describe below. 
We have 

V($) cx 2A cdt $- i$ 3 , (28) 
and the cubic potential is shown in Fig. 3. It has a local minimum £q determined by 

V'(£ ) = =>■ l = -V2A cd t. (29) 

Thus the "classical" limit of the matrix integral with potential V($>), when only the 
minimum plays a role, leads to the following expectation value: 

4 (tr 1 \ = ^== = lim W cdt (Z cdt , A cdt , G,), (30) 

so in this way one can use this matrix model very explicitly to obtain the classical 
limit. As we will show in the next Section even this classical limit has a non-trivial 
representation as a sum over certain random surfaces. 

Geometric interpretation 

Let us define some geometric objects, related by Laplace transformations, W\,g. (£) 
and W\ g s (x). We have now a little abuse of notation. Above we used capital letters 
for continuum, dimensionful variables and small letters for dimensionless (lattice) vari- 
ables. Now we will use A and g s x for continuum, dimensionful variables. In the end 
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Figure 4: The graphic representation of the integral equation satisfied by the full disk 
amplitude, eq. (32). 



they will be identified with the continuum variables A cc jt, G s and Z c d t etc., in the same 
way as the object W\,g s {%) in the end will be identified with W c dt(Z c dt, A cc jt, G s ). 
However, we do it in order to stress that we are now starting from scratch with a geo- 
metric theory which in principle has nothing to do with the matrix model above. And 
in fact that was how these "quantum geometric" objects to first found and analyzed 



The objects are intended to represent the disk amplitude in the theory of two-dimensional 
quantum gravity based on causal dynamical triangulated random surfaces (CDT). The 
idea is to start be summing over all random surfaces with a boundary of length I which 
admit a proper time foliation. The action used is just the area action (i.e. the cosmolog- 
ical term, since the Einstein term is trivial in 2d as long as we do not allow topology 
change. And we will not allow that presently. Thus the topology of the surface is just 
the trivial topology of the disk.) We denote this sum the disk amplitude w{ (£). It is 
called the CDT disk amplitude. If we add a boundary cosmological constant x at the 
boundary, we should add a boundary action x £, and thus the Laplace transformation 
(31) can be viewed as changing the path integral over surfaces with a fixed boundary 
length to a path integral where we also integrate over all boundary lengths and instead 
keep fixed a boundary cosmological constant x. 

We now allow a larger class of surfaces by allowing branching, i.e. we allow the 
spatial surface at a proper time t to split in two, the splitting assigned a weight g s . g s is 
clearly like a string coupling constant. The process is shown in Fig. 4. We are allowing 
the splitting of a spatial universe in two, but presently we do not allow for topology 
change of the two-dimensional surface, so we do not allow the two universes to join 
again, since that would create a handle and then change the two-dimensional topology. 
Had we allowed it, the total coupling constant associated with this process would have 
been g^, one factor for splitting, one factor for joining, like in string theory. The un- 
shaded disk amplitude is Wjf* (£), while the full disk amplitude is denoted W\ t g s (£) 
and is shown as the shaded graph on the lhs of the equality sign. 



[12]-[13]. 




(31) 
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Figure 5: The graphic illustration of the geometries which contribute to the propagator 
G\,g s (x,y;t) ineq. (34). 



Fig. 4 is a graphic representation of the following integral equation: 



A 

OO OO 

g s J dt 
o o 



(32) 



d£ 2 (ii+i2)Gf ] (xJ 1 +e 2 ;t)W x ^(ii)W^ gs (£ 2 ) 



In this equation the object G^ (£i ,l2]t) denotes the "propagator" in CDT. It is defined 
in analogy with W> (£). We sum over all two-dimensional geometries where we have 
a spatial entrance loop of length l\ and a spatial exit loop of length 1%, with the further 
constraint that all point on the exit loop is separated a geodesic distance t from the 
entrance loop. Again we assume that all points separated a geodesic distance t' < t 
from the entrance loop form a connected one-dimensional space. G[ a) {x,£';t) denotes 
the Laplace transform of G^ (£, £'; t) with respect to I. 

In the same way as we generalized the geometries which entered into path integral 
defining {£), by allowing space to split, and in this was obtained W\ gs {£), we can 
allow for the spatial hyper-surface at any time t' < t to separate in two. One of these 
will then be connected to the exit loop while the other will eventually vanish in the vac- 
uum. This is shown in Fig. 5. We denote this generalized propagator G\, 3s (£i,£ 2 ; t)- 
For both G x , 9B {h,h;t) and G x 0) (£ 1 ,£ 2 ;t) we introduce, again in analogy with the 
definitions for W\, Bs (£) and Wj®(£), the Laplace transforms: 



G\,g s {x,y;t) 



d£i 



M 2 e" 



(33) 



as well as the hybrid forms G\ i9s (x, £ 2 ;t) and G\ i9s (£i, y; t). Here x and y denotes 
boundary cosmological constants at the entry and exit boundaries. 

The shaded parts of graphs in Fig. 5 represent the full, ^-dependent propagator and 
the full g s -dependent disc amplitude, and the non-shaded parts the CDT propagator 
where g s = 0. In all four graphs, the geodesic distance from the final to the initial loop 
is given by t, while the "baby-universes" which end in the vacuum can terminate their 
life at any positive time after they have been created, also at a time larger than t. 
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From Fig. 5 one can write down an integral equation much like eq. (32) for W\ j9e (x). 
However it is convenient the differentiate this equation with respect to t and we then 
obtain: 



a£ Ql G ^s V\ t) = -^ [(a(x 2 - A) + 2g s a s aP^W^g, (a)) G x , 9s (x, y; t)] , 

(34) 

In this equation we have explicitly assumed that we have some kind of regularized the- 
ory, on a lattice, say. The equation is then first written in terms of the dimensionless 
lattice variables and the translated into continuum notation by inserting the relation 
between the dimensionless variables and their continuum counterparts, in this way in- 
troducing the lattice cut-off a. This cut-off will then appear with a power determined 
by the dimension of the continuum variables, except for a subtlety related to W, to be 
explained now. The lattice cut-off a is assigned length dimension 1, the time t assigned 
the unspecified length dimension e, and the disk amplitude W the length dimension 
—77. We assume x, as the coupling constant conjugate to the length £ has length dimen- 
sion — 1. We leave the length dimension of g s unspecified as —6. 

Denote the dimensionless lattice variables and observables by t reg , W reg etc. We 
thus have t reg = t/a e . Similarly with W and W reg : W reg — Wa v , except that 
we allow for the possibility that W reg is not scaling when a — >• (i.e. that 77 > 0). 
Although at first sight a little strange this was precisely what happened in the ordinary 
scaling limit of the matrix models, as we noted in the discussion above (see eq. (11) 
which shows that in the ordinary scaling limit we have 77 = 3/2). Thus we will allow 
for this possibility. We can summarize as follows: 

W reg ► aPW x {x), 77 <0, (35) 

a— >0 

treg > t / fl£ , £ = 1. (36) 



Wreg > const. + a? W x (x), 77 = 3/2 (37) 

a—>0 

Ueg ► t/cf, E = 1/2, (38) 

a— ^0 

The values of the exponents e and 7/ are written to the right in equations (35)-(38). 
They tell us that when 7/ < 0, i.e. when W\ t9s (x) scales, then s = 1, while if 7/ > 
then £ = 1/2 and r\ = 3/2. The reason these exponent are uniquely determined is that 
we have the geometric picture shown in Fig. 6, which couples W and G and leads to 
consistency relations for the scaling of W: 



dW±gM 



dt dlGx, g .{x,t;t)Wx, g .(t). (39) 



We know the scaling dimension of G\. gs (x, £; t). It is zero. It all goes back to the fact 
that G as a propagator has to satisfy the composition rule: 

/•oo 

G(e 1 ,£ 2 ;t 1 +t 2 )= / d£G(£x,i;h)G(£,£ 2 ;h), (40) 
Jo 
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Figure 6: The differentiation of the disk amplitude with respect to the cosmological 
constant leads to the disk amplitude with one marked point. Such a geometry has the 
unique decomposition shown in the figure, where the loop has a geodesic distance t 
from the boundary loop. The equation graphically represented on the figure is (39). 

valid for both G x ,g s (4,4; *) and G^ 0) (4, 4; t). This means that (7(4,4; t) has to 
scale like a^ 1 and the Laplace transform G(x,£;t) thus as a°. Combined with that 
fact that the cosmological constant A has length dimension —2, we are led to right 
side values of e and r) in eqs. (35)-(38). It is a beautiful example of the constraints 
imposed by quantum geometry, and it is remarkable that one is able to derive the non- 
trivial, ordinary matrix model value rj = 3/2 in eq. (11) from such simple geometric 
considerations. 

If we choose the solution (35)-(36) then we can solve the geometric equations (32) 
and (34). A glance on eq. (34) shows that //the term involving g s is going to play a 
role we have to take S = 3, the same result as in the new scaling limit of the matrix 
models. If we make this choice the solution for W\ t9s (x) is precisely the one given 
by the new scaling limit of the matrix models. Thus the new scaling limit has indeed 
a geometric representation in terms of random surfaces, and even a regularized lattice 
representation where a lattice spacing a is taken to zero. All this is discussed in detail 
in [12]-[13]. The main characteristic is here that the geodesic distance (the proper time 
t) has canonical scaling dimension, identical to that of space and there is a smooth limit 
g s — > where one obtains the original CDT solution [9]. 

If we choose the solution (37)-(38) and solve the geometric equations (32) and (34) 
we find a completely different solution. It is characterized by a different scaling of 
the geodesic distance or proper time t. We see that eq. (34) is only non-trivial if we 
choose the dimension of g s to be zero, and the first term on the rhs of eq. (34) is then 
irrelevant. The equation is thus reflecting an excessive branching off of baby universes. 
There is nothing else! Thus, looking at Fig. 4 and 5, we have no unshaded parts of the 
graphs. They simply play no role in the scaling limit where the lattice spacing a 0. 
The typical geometry which appears in the path integral will be very fractal and will 
have Hausdorff dimension dh = 4, not dh = 2 as one would expect from a "nice" two- 
dimensional geometry. This is described in detail in [15]. The wild branching of baby 
universes is illustrated in Fig. 7. This limit with dh = 4 corresponds to the "ordinary" 
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Figure 7: Part of a "typical" triangulation of the disk which one will encounter in 
the "ordinary scaling" limit of the matrix models, where equations (37) and (38) are 
satisfied. The dashed and dotted lines represent two "spatial" curves separated by one 
(proper)-time step, plus all the baby universes which are cut off at the this time step. 

scaling limit of the matrix model. 

This is in sharp contrast to the solution provided by conditions (35)-(36). There the 
typical geometry has dh = 2 and the total number of baby universes created is finite. 
This implies that a typical geometry will look as shown in Fig. 8. In Fig. 8 we have 
drawn an amplitude which is more complicated than the disk amplitude (we have two 
entrance loops, and the surface also have a handle, i.e. it is a higher genus surface. This 
gives us the opportunity tO emphasize that while we have mainly discussed geometries 
of the simplest topology, the whole discussion of the new scaling generalizes to any 
genus surface and with any number of boundary loops, and it matches precisely the 
1 /N 2 expansion in the matrix models of expectation values of multiple trace operators 
[12]-[14]: 

w(zi, ...,Z n )= (— tr 1 —— ■ ■ ■ — tr connected- (41) 

N z\ — <p Jy z n — <p 

These multiple trace operators are obvious generalizations of the single trace operator 
defined in eq. (7). 
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Figure 8: A "typical" geometry (in the continuum limit) of the disk which one will 
encounter in the "new" scaling limit of the matrix models, where equations (35) and 
(36) are satisfied. We have show a surface where the topology of the surface is also 
changed from that of a simple disk to a disk with a handle, see discussion in the main 



Unfinished stuff 

For the ordinary matrix models we have a description of conformal matter coupled to 
2d quantum gravity by multicritical one-matrix models and by two-matrix models. 

Similarly it is easy to couple matter to the "plain" CDT model. Ising models and 
multiple Pott models coupled to CDT have been studied numerically [16]. 

Now that we have a matrix model description of the generalized CDT models it 
seems natural to apply the same technique as was applied for the ordinary matrix mod- 
els and in this way use matrix models to study the matter systems defined on the CDT- 
like set of random surfaces. From the computer simulations referred to, it seems that 
one obtain the fiat space-time exponents. It would be very interesting if one could 
obtain a simple proof of the any conformal field theory exponent from a matrix inte- 
gral. It would provide us with an explicit realization of these critical systems, even at a 
regularized level. Work in this direction is in progress. 
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